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Remarks on the discretization of physical momenta in lattice QCD* 
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The calculation on the lattice of cross-sections, form-factors and decay rates associated to phenomenologically 
relevant physical processes is complicated by the spatial momenta quantization rule arising from the introduction 
of limited box sizes in numerical simulations. A method to overcome this problem, based on the adoption of two 
distinct boundary conditions for two fermions species on a finite lattice, is here discussed and numerical results 
supporting the physical significance of this procedure are shown. 



1. The problem 

In the formulation of quantum field theories on 
a lattice the introduction of a finite volume is un- 
avoidable when numerical simulations are used as 
tool of investigation. As a consequence of the lim- 
ited box size, spatial momenta come out to be 
quantized accordingly to the choice of the bound- 
ary conditions. 

The momentum quantization represents a se- 
vere limitation in various phenomenological ap- 
plications by making impossible, in the great ma- 
jority of the cases, to calculate cross sections, 
decay-rates and form-factors in the interesting 
kinematical region or at the physical values of the 
masses of the particles involved in the processes. 
In this talk it is discussed a method, previously in- 
troduced in ref. PP , that allows to overcome these 
difficulties. The idea consists in using the depen- 
dence of the momentum quantization condition 
upon the choice of the boundary conditions (BC) 
and in requiring different fermion species to sat- 
isfy different BC. 

2. The way out 

In order to clarify the dependence of the mo- 
mentum quantization rule upon the choice of the 
boundary conditions let us first consider the case 
of a particle satisfying periodic boundary condi- 
tions (PBC). For a fermionic field ip(x) on a 4- 
dimensional finite volume of topology T x L 3 with 



PBC in the spatial directions one has 

ip(x + Bi L) = ip(x) , i = 1,2,3 (1) 

This condition can be re-expressed by Fourier 
transforming both members of the previous equa- 
tion 
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rP(p) = / d A p e"** 4,(p) (2) 



and implies 

e^ L = 1 = 
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where the n^'s are integer numbers. The authors 
of [213141 5 6 7 8 9 have considered a generalized 
set of boundary conditions, that here we call 9- 
boundary conditions (f?-BC), depending upon the 
choice of a topological 3-vector 6 



ip(x + ei L) = e Wl ip(x) 



1,2,3 



(4) 



The modification of the boundary conditions af- 
fects the zero of the momentum quantization rule. 
Indeed, by re-expressing equation 10} in Fourier 
space, as already done in the case of PBC in equa- 
tion ©, one has 



*(P»-T") L = 1 
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Pi= L 



(5) 
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It comes out that the spatial momenta are still 
quantized as for PBC but shifted by an arbi- 
trary continuous amount (Qi/L). The physical 
significance of this continuous shift in the allowed 
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momenta becomes manifest, from the theoretical 
point of view, by realizing that eq. Q is the con- 
dition satisfied by the electron wavefunctions in 
a periodic solid crystal accordingly to the well 
known Bloch's theorem. Indeed, quarks on a lat- 
tice experience a periodic potential (the gauge 
fields) exactly as the electrons do in a crystal and 
the results of the Bloch's theorem apply straight 
forwardly also to the lattice fields. 

The generalized ^-dependent boundary condi- 
tions of equation J3J can be implemented by mak- 
ing a unitary Abelian transformation on the fields 
satisfying 0-BC 



ijj(x) 



U (9,x)ip(x) = e- 1 ^ ip(x) (6) 



As a consequence of this transformation the re- 
sulting field satisfies periodic boundary condi- 
tions but obeys a modified Dirac equation 

J2$( x ) D e{x,v) (7) 

where the ^-dependent lattice Dirac operator 
Dg(x,y) is obtained by starting from the pre- 
ferred discretization of the Dirac operator and by 
modifying the definition of the covariant lattice 
derivatives (see ref. pQ for details). 

3. Numerical tests 

In order to prove numerically that the term 
9/L acts as a true physical momentum, one can 
study the energy of a meson made up by two 
different quarks with different 9-BC for the two 
flavors. In the following we work in the O(a)- 
improved Wilson— Dirac lattice formulation of the 
QCD within the Schrodinger Functional formal- 
ism but, we want to stress that the use of 9-BC 
in the spatial directions is completely decoupled 
from the choice of time boundary conditions and 
can be profitably used outside the Schrodinger 
Functional formalism, for example in the case of 
standard periodic time boundary conditions. Let 
us consider the following correlators 

f$(9;x ) = 
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Figure 1. Effective energies E^^(0, a; xq) at fixed cut- 
off. The resuits correspond to the simulation done at 
P = 6.211 with r mf G/ = 0.655 and r m^ GI = 0.354. 
Similar figures could have been shown for other combi- 
nations of the simulated quark masses and for the other 
values of the bare coupling. 
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where i and j are flavor indices, all the fields sat- 
isfy periodic boundary conditions and the two fla- 
vors obey different ^-modified Dirac equations, as 
explained in eq. Q. In practice it is adequate to 
choose the flavor i with 9 = 0, i.e. with ordi- 
nary PBC, and the flavor j with 9^0. After the 
Wick contractions the pseudoscalar correlator of 
equation JSJ reads 

= Y E Tr ^-(0; z, x) 75 5 l (0; x, y)) (9) 

y,z,x 

where S(9;x,y) and S(0; x, y) are the inverse of 
the ^-modified and of the standard Wilson-Dirac 
operators respectively. Note that the projection 
on the momentum 9/L of one of the quark legs 
in equation it is not realized by summing on 
the lattice points with an exponential factor but 
it is encoded in the ^-dependence of the modified 
Wilson-Dirac operator and, consequently, of its 
inverse S(9;x,y). This correlation is expected to 
decay exponentially at large times as 
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where, a part from corrections proportional to the 
square of the lattice spacing, is the physical 
energy of the mesonic state. After the continuum 
extrapolations one has to recover the expected 
relativistic dispersion relations 

4 = ^ + (z) (") 

where My is the mass of the pseudoscalar meson 
made of a i and a j quark anti-quark pair. 

All the numerical results are obtained in the 
quenched approximation of the QCD. We have 
done simulations on a physical volume of topology 
T x L 3 with T = 2L and linear extension L = 
3.2 ro, where ro is defined in [Hi]- All the missing 
parameters of the simulations are given in table 1 
of ref. [J. 

Setting the lattice scale by using the physical 
value r = 0.5 fm, the expected values of the 
physical momenta associated with our choices of 9 
are calculated according to the following relation 

! 0.000 
° Q f 3 l GeV (12) 
0.650 

where L ~ 1.6 fm. These values have to be com- 
pared with the value of the lowest physical mo- 
mentum allowed on this finite volume in the case 
of periodic boundary conditions, i.e. \p\ ~ 0.785 
GeV. 

At fixed cut-off, for each combination of fla- 
vor indices and for each value of 9 we have ex- 
tracted the effective energy from the correlations 
of eq. IJSJ, fp(8;x ). In fig. ^ we show this 
quantity for the simulation performed at (3 — 
6.211 corresponding to ro m^ GI = 0.655 and 
ro rri2 GI = 0.354, for each simulated value of 9. 
As can be seen the correlations with higher val- 
ues of If? | are always greater than the correspond- 
ing ones with lower values of the physical mo- 
mentum, a feature that will be confirmed in the 
continuum limit. Being interested in the ground 
state contribution to the correlation of eq. ©, we 
have averaged the effective energies over a ground 
state plateau of physical length depending upon 
the quark flavors and we have extrapolated the 
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Figure 2. Continuum dispersion relations. The data cor- 
respond to different combinations of the simulated quark 
masses and reproduce very well the expected theoretical 
behavior, i.e. straight lines having as intercepts the meson 
masses and as angular coefficients one (see eg. 1111 , 

resulting quantities, E lJ (9,a), to the continuum 
(see ref. pQ for details). 

The continuum results verify very well the dis- 
persion relations of equation l|ll|) as can be clearly 
seen from fig. [3 in which the square of E^(9) for 
various combinations of the flavor indices is plot- 
ted versus the square of the physical momenta 
\9\/L. The plotted lines have not been fitted but 
have been obtained by using as intercepts the sim- 
ulated meson masses and by fixing their angular 
coefficients to one. 

REFERENCES 

1. G. M. de Divitiis, R. Petronzio and N. Tantalo, 
Phys. Lett. B595 (2004) 408, hep-lat/0405002 

2. K. Jansen et al., Phys. Lett. B372 (1996) 275, 
|hep-lat/9512009| 

3. A. Bucarelli et a l, Nucl. Phys. B552 (1999) 379, 
|hep-lat/ 9808005 

4. Zeuthen-Rome / ZeRo, M. Guagnelli et al., Nucl. 
Phys. B664 (2003) 276, hep-lat/0303012 

5. PF. Bedaque, (2004), nucl-th/0402051 

6. D.J. Gross and Y. Kitazawa, Nucl. Phys. B206 
(1982) 440. 

7. J. Kiskis, R. Narayanan and H. Neuberger, Phys. 
Rev. D66 (2002) 025019, hep-lat/0203005 

8. J. Kiskis, R. Narayanan and H. Neuberger, Phys. 
Lett. B574 (2003) 65, hep-lat/0308033 

9. A. Roberge and N. Weiss, Nucl. Phys. B275 
(1986) 734. 

10. R. Sommer, Nucl. Phys. B411 (1994) 839, 
|hep-lat/9310022| 



